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silicon lenses for millimeter wavelengths 


R. Datta, 1 * C. D. Munson , 1 M. D. Niemack , 23 J. J. McMahon , 1 J. Britton , 3 E. J. Wollack , 4 
J. Beall , 3 M. J. Devlin , 5 J. Fowler , 3 P. Gallardo , 6 J. Hubmayr , 3 K. Irwin , 3 L. Newburgh , 7 
J. P. Nibarger , 8 L. Page , 7 M. A. Quijada , 4 B. L. Schmitt , 5 S. T. Staggs , 7 

R. Thornton , 9 and L. Zhang 7 

University of Michigan Physics Department, 450 Church St., Ann Arbor, Michigan 48109, USA 

2 Cornell University Physics Department, 109 Clark Hall, Ithaca, New York 14853, USA 

3 National Institute of Standards and Technology, 325 Broadway, Boulder, Colorado 80305, USA 

4 NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, Maryland 20771, USA 

5 University of Pennsylvania Department of Physics and Astronomy, 209 South 33rd St., Philadelphia, Pennsylvania 19104, USA 

6 Departamento de Astronomia y Astrofisica, Facultad de Fisica, Pontificia Universidad Catolica, Casilla 306, Santiago 22, Chile 

7 Princeton University Department of Physics, Jadwin Hall, P.O. Box 708, Princeton, New Jersey 08544, USA 

8 Boulder Micro-Fabrication Facility, National Institute of Standards and Technology, 325 Broadway, 

MS 817.03, Boulder, Colorado 80305, USA 

9 West Chester University of Pennsylvania Department of Physics, 700 South High St., West Chester, Pennsylvania 19383, USA 

"Corresponding author: dattar@umich.edu 

Received 26 July 2013; revised 24 October 2013; accepted 25 October 2013; 
posted 25 October 2013 (Doc. ID 194154); published 16 December 2013 


The increasing scale of cryogenic detector arrays for submillimeter and millimeter wavelength astrophys- 
ics has led to the need for large aperture, high index of refraction, low loss, cryogenic refracting optics. 
Silicon with n = 3.4, low loss, and high thermal conductivity is a nearly optimal material for these 
purposes but requires an antireflection (AR) coating with broad bandwidth, low loss, low reflectance, 
and a matched coefficient of thermal expansion. We present an AR coating for curved silicon optics com- 
prised of subwavelength features cut into the lens surface with a custom three-axis silicon dicing saw. 
These features constitute a metamaterial that behaves as a simple dielectric coating. We have fabricated 
silicon lenses as large as 33.4 cm in diameter with micromachined layers optimized for use between 125 
and 165 GHz. Our design reduces average reflections to a few tenths of a percent for angles of incidence 
up to 30° with low cross polarization. We describe the design, tolerance, manufacture, and measurements 
of these coatings and present measurements of the optical properties of silicon at millimeter wavelengths 
at cryogenic and room temperatures. This coating and lens fabrication approach is applicable from cen- 
timeter to submillimeter wavelengths and can be used to fabricate coatings with greater than octave 
bandwidth. © 2013 Optical Society of America 

OCIS codes: (000.2190) Experimental physics; (220.0220) Optical design and fabrication; (220.3630) 

Lenses; (310.1210) Antireflection coatings; (160.1245) Artificially engineered materials; (310.6628) Sub- 
wavelength structures, nanostructures. 
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1. Introduction 

The development of large format superconducting 
detector arrays for millimeter and submillimeter 
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astrophysics (e.g., [1—3] ) has driven the need for high- 
throughput optical designs that maintain diffraction 
limited performance across the arrays (see [4] for a 
review). Silicon is an excellent material for these ap- 
plications due to its high index of refraction (n — 3.4), 
low loss-tangent (tan <5 < 7 x 10" 5 at 4 K), and rela- 
tively high thermal conductivity ( k ~ 200 Wm' 1 K -1 
at 4 K [5] , which is only a factor of 50 lower than high 
quality oxygen free high conductivity (OFHC) copper 
[6] and orders of magnitude higher than plastic). The 
essential development required to realize optical de- 
signs using silicon optics at millimeter wavelengths 
is an appropriate antireflection (Aft) coating to 
mitigate the 30% (-5 dB) reflective loss from each 
optical surface. 

AR coatings consist of one or more dielectric layers 
placed on the surface of refractive optical elements. 
The indices of refraction and thicknesses of the AR 
layers are chosen such that the reflections from the 
vacuum-AR interface and from the AR-substrate 
interface interfere and cancel. If the optical path 
length through a single-layer AR coating is one quar- 
ter wavelength, the refractive index of the coating is 
s/n (the geometric mean of the substrate index and 
free space), and the coating is free from dielectric 
losses, the cancellation is perfect at one particular 
frequency. Applied to silicon, such quarter-wave coat- 
ings reduce reflections below 1% over 1.25 : 1 band- 
width at normal incidence. 

Wider bandwidths and larger ranges of angles of 
incidence can be accommodated by adding additional 
layers to the AR coating to form a multilayer coating. 
These coatings require layers with a number of 
different refractive indices. Cryogenic applications 
require that the coefficients of thermal expansion of 
the AR layers are sufficiently matched to that of the 
lens substrate to prevent damage upon cooling. The 
dielectric loss and birefringence of the AR coating 
material must also be carefully controlled. Locating 
materials with these properties poses a significant 
challenge to implementing cryogenic wide band AR 
coatings for silicon. One avenue to solve this problem 
is to engineer materials with the required dielectric 
constant by cutting subwavelength features into the 
lens surface [7,8] . Engineered optical materials with 
properties determined by their detailed geometric 
shape, size, and orientation are referred to as “arti- 
ficial dielectric metamaterials.” 

These metamaterial coatings consist of layers of 
subwavelength holes [9], posts, or grooves [ 10 ] cut 
into the substrate to be coated. This approach allows 
the index of the AR layer to be precisely tuned by 
adjusting the geometry of the machined features, 
naturally solves the thermal expansion problem, and 
has loss lower than that of the substrate. At cryo- 
genic temperatures silicon has a dielectric loss tan- 
gent ~100 times lower than typical plastics like 
Cirlex, which have previously found use as single- 
layer AR coatings [ 11 ] . Rosen et al. [12] developed an 
epoxy-based multilayer broadband AR coating for 
millimeter-wave cryogenic optics. The advantages 


of silicon metamaterial AR layers include precisely 
controlled indices of refraction, inherently matched 
thermal expansion, and significantly reduced dielec- 
tric losses. 

This approach has been successfully applied to plas- 
tic millimeter-wave optics through direct machining 
and silicon optics for visible wavelengths through 
patterning of geometric structures which behave as 
an AR coating similar to that found on a moth’s eye 
[13—15] . At millimeter wavelengths Zhang et al. [ 16 ] 
fabricated an alternative class of artificial dielectric 
AR coatings for silicon composed of metal resonant 
structures supported by a plastic substrate that had 
higher losses than the AR layers described here. 
Schuster et al. [ 17 ] fabricated single-layer silicon mi- 
cromachined artificial dielectrics, and Han et al. [ 18 ] 
have used silicon immersion grating technology to 
fabricate planar artificial dielectrics for terahertz ra- 
diation. Our work is the first demonstration of broad- 
bandwidth metamaterial AR layer on silicon optics 
with finite curvature for millimeter wavelengths. 

The AR coated lenses presented here have been 
developed for the ACTPol [ 19 ] project, a polarization 
sensitive receiver for the Atacama Cosmology Tele- 
scope (ACT, [20]). In Section 2 we present the design 
requirements for the ACTPol optics as they apply to 
this work. The design of the AR layers is described in 
Section 3. In Section 4 we describe the ACTPol lens 
design. Cryogenic measurements of the dielectric 
properties of silicon samples taken from 45 cm diam- 
eter ingots at millimeter wavelengths are presented 
in Section 5. In Section 6 we describe the fabrication 
process of the AR coatings on the lenses and present 
reflection measurements of a finished lens. We con- 
clude with a brief discussion of the range of appli- 
cability of this coating technique. 

2. Requirements 

The performance requirements for the AR coated 
silicon lenses described in this work are dictated 
by the scientific goals of the ACTPol project [ 19 ] . This 
instrument is designed to detect the faint polarized 
signals from the cosmic microwave background 
(CMB) using three detector arrays fed by indepen- 
dent reimaging optics. Two of the arrays are horn 
coupled polarimeters for the 150 GHz band (pass- 
band from 125 to 165 GHz), while the third array will 
be a multichroic array including both 150 and 
90 GHz bands (“90 GHz” passband from 80 to 
110 GHz) [21]. Extensive design studies demon- 
strated high index of refraction (n > 3) lenses were 
needed to achieve high optical quality across the re- 
quired field of view. 

In this paper we focus on the AR coatings for the 
150 GHz band lenses. These are conceptually similar 
to those which will be used for the broader-band mul- 
tichroic lenses. The bandwidth requirement for the 
150 GHz band coatings is to have >99% transmission 
between 125 and 165 GHz. The left panel of Fig. 1 
shows the reimaging optics for a 150 GHz array. The 
optics include three plano-convex silicon lenses with 


8748 


APPLIED OPTICS / Vol. 52, No. 36 / 20 December 2013 


lens reimaging optics 


distribution of angles of incidence 



angle of incidence (deg) 



Fig. 1. Left: Ray diagram of one set of ACTPol reimaging optics, 
which includes three silicon lenses feeding a detector array. Right: 
Histogram of the angles of incidence of rays at the surfaces of these 
three silicon lenses. 


diameters up to 33.4 cm. The figure shows that rays 
passing through the optical system refract over a 
wide range of angles of incidence. This is quantified 
in the right-hand panel of Fig. 1, which shows that 
the distribution of angles of incidence is centered 
near 17° and that more than 96% of the rays have 
angles of incidence <30°. This sets the requirement 
that the coatings must be optimized to minimize re- 
flections for angles of incidence between 0° and 30°. 

Since ACTPol is a polarization sensitive experi- 
ment, low cross polarization is another requirement. 
Studies of polarization systematics (e.g., [22]) sug- 
gest that the CMB temperature to polarization 
leakage must be controlled to better than 1%, which 
corresponds to a requirement that differences in the 
transmission for the two polarizations be 0.5%. 

To reduce the thermal emission, all optics are 
cooled to 4 K or below. Thus the AR coating must 
be able to withstand cryogenic cycling. 

As described in Section 5, high resistivity silicon, 
which is available in boules up to 45 cm in diameter, 
offers a cryogenic loss tangent tan 5 <7 x 10 -5 , an 
index of refraction of n — 3.4, and a relatively high 
thermal conductivity, making it an ideal material 
for this optical design. Our approach of directly 
machining metamaterial AR layers into the lens sur- 
faces guarantees that the layers have low dielectric 
losses and coefficients of thermal expansion that are 
inherently matched to that of the silicon lenses. 

3. AR Coating Design 

The AR coating is comprised of two metamaterial 
layers. Each layer consists of an array of square 
pillars on a square grid cut using silicon dicing 
saw blades by making two sets of parallel trenches 
rotated 90° relative to each other. It has been shown 
that any electromagnetic structure with an axis of 
n-fold (n > 2) rotational symmetry must have polari- 
zation independent reflectance and zero cross polari- 
zation when light is incident along this axis [23], 
Hence the rotational symmetry of these layers leads 
to zero cross polarization at normal incidence and 
low birefringence at oblique incident angles. Figure 2 
shows a cross section and isometric view of this 
geometry which is parameterized by the depths d\, 
d 2 and kerf widths k\, k 2 of the inner and outer layer, 
respectively, and a pitch p. The outer layer is cut 
by using a blade of width k 2 and the inner layer is 


cross-section 




Fig. 2. Geometry of the two-layer metamaterial AR coating. The 
left side shows a cross section through the center of the pillars 
indicating the design parameters discussed in the text. The right 
side shows an isometric view of the structure resulting from 
making cuts along two perpendicular directions. 


cut by making a second set of cuts to greater depth 
with a thinner blade of width k\. We introduce the 
volume fill factor Vf — (1 - k/p) 2 to facilitate com- 
parison to published analytic calculations based on 
the second order effective medium theory developed 
by Rytov [ 24 ] for the effective dielectric constant in 
the quasi-static limit (p <c X) [15,25-27] . 

The metamaterial layers can be treated as a vol- 
ume distribution of small electromagnetic scatterers 
characterized by electric and magnetic polarizability 
densities [7] . When the pitch is small compared to 
the shortest wavelength of interest, the fields in the 
layer are homogenous [8] and one can define an 
effective dielectric constant which can be used to pa- 
rameterize the propagation properties of the media. 
This effective dielectric function for the layer, e® ff , is a 
function of the density of scatterers in the layer 
[28,29], parametrized by the volume fill factor. As 
the feature size of a composite medium approaches 
a significant fraction of the radiation wavelength, 
the effective dielectric function becomes frequency 
dependent [ 30 ] . As the wavelength is further reduced 
the artificial dielectric structure no longer appears 
homogenous and these simple quasi-static consider- 
ations must be augmented to adequately model its 
behavior. In going beyond the quasi-static limit 
resonant effects, diffraction and scattering can occur 
in the artificial dielectric structure [9,31—33] . The AR 
coatings described here, which operate between the 
quasi-static and diffractive regimes, rely on numeri- 
cal simulations to estimate the relation between the 
volume fill factor and the effective refractive index 
of the metamaterial layer. 

The values of d 1; d 2 , k x , k 2 , and p were chosen us- 
ing a three-step design process. A preliminary design 
was carried out using a classical analytic model 
consisting of sheets of dielectric material each with 
a constant thickness and index of refraction [34], 
These were varied to minimize the reflectance across 
the band yielding targets for the electrical thickness 
and index of refraction of our two-layer coating. 

The second step was to translate the index and 
thickness from the simple model into the pitch, kerfs, 
and depths for our geometry. This step required 
knowledge of the relation between the effective index 
of refraction of an array of square pillars and the 
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Fig. 3. Comparison of the reflectance calculated from an HFSS 
numerical simulation of a pillar geometry to fits based on a model 
of a simple dielectric layer. The geometry is specified in the inset in 
the top. In the simulated geometry, p = 400 pm, £ phys = 220 pm, 
k = 40 pm, and re si = 3.38. A comparison between fits where only 
the index n is free and where both the index n and thickness t are 
varied is shown. The HFSS simulation and the best fit curve differ 
by 5% at «220 GHz (1.36 mm) for the analytic fit varying n and t 
and at «180 GHz (1.67 mm) for the analytic fit varying only n. This 
corresponds to the minimum wavelength, 1 5% , (and a correspond- 
ing maximum frequency f 5% ) for the specified pitch, p, below 
(above) which the layer no longer behaves as a simple dielectric 
layer. 

pitch and kerf of that layer. These relations were de- 
termined by fitting analytic models to high frequency 
structure simulator (HFSS, [35]) simulations of single 
layers with a wide range of volume fill factors for a few 
different choices of pitch. The HFSS simulations were 
carried out using Floquet ports and master-slave 
boundary conditions to model a box containing a sin- 
gle post of the metamaterial as an infinite periodic 
two-dimensional array of features [ 35 ] . 

Figure 3 shows a representative simulation of a 
single-layer pillar geometry with insets showing the 
simulated geometry and the simple dielectric layer 
modeled in the fit. The fit constrained the effective 
index n eS , the effective electrical thickness t eS , and 
the maximum frequency to which our metamaterial 
layer behaves as a simple dielectric. For the purposes 
of our fit we define this maximum frequency as the 


point where the simulated reflectance and the ana- 
lytic fit disagree at more than 5% absolute reflectance. 
We refer to this frequency as f 5% and the correspond- 
ing wavelength A$%. This indicates the transition be- 
tween specular and diffractive behavior of the layer. 

A second fit which fixed the electrical thickness to 
the physical thickness of the pillars is also shown in 
Fig. 3. This fit implies a lower f 5% . The discrepancy 
between the best fit thickness t e ff and the physical 
thickness t phys of the metamaterial layer arises since 
the geometric structure of the pillars results in fring- 
ing of fields at the junction between layers, causing a 
perturbative shift between the position of the physical 
interface and the effective location of the junction’s 
electrical reference plane. An analogous effect is en- 
countered in the design of metallic waveguide struc- 
tures [32], which can be analytically treated. 

Figure 4 shows the results from fitting a number of 
HFSS simulations with different pitch and fill factors 
to the analytic model. Plots include the breakdown 
wavelength A 5% , the effective index n e ff, and the elec- 
trical thickness of a layer with index ,/ftsi expressed 
as a ratio of the thickness of a quarter-wave homog- 
enous dielectric layer and a quarter-wave metamate- 
rial layer. 

The recovered A 5% is in good agreement with the 
analytic condition [9,33] to prevent diffraction from 
a grating array, 


P < 


A 

(n s + n, sin e t ) ' 


( 1 ) 


where p is the pitch, A is the wavelength, n s is the in- 
dex of the substrate, i%i is the index of the incident 
medium, and (9, is the angle of incidence. This analytic 
expression shows that the pitch must be chosen based 
on the indices of the adjacent media and must be 
smaller than implied by the 1 5% calculated at normal 
incidence to minimize diffraction at oblique angles of 
incidence. Consideration of manufacturing cost and 
mechanical robustness of the two-layer AR design fa- 
vors the largest possible pitch, as this choice minimizes 
the number of cuts required to cover a given area and 




Fig. 4. Summary of results from HFSS simulations of single layers of pillars with varying geometry. Left: Minimum wavelength for 
which the single pillar layer is fit by an analytic model to 5% accuracy as a function of pitch p at normal incidence. Center: effective index of 
refraction ra eff as a function of the volume fill factor Vf obtained by fitting the simulated reflectances up to the minimum wavelength 1 5% = 
1360 pm for a pitch p = 400 pm, 40 pm. The line labeled “Biber model” is a prediction from the analytic Biber model [25], which corre- 
sponds to the quasi-static circuit approximation for the media. Right: Thickness t e g for a slab dielectric quarter-wave layer divided by the 
thickness i pllys of pillars forming an effective quarter-wave coating as a function of p /X. The index of the simulated dielectric layer was held 
fixed at y/n^ for these simulations. A quadratic fit is plotted to guide the eye. 
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makes the pillars larger and therefore stronger. A 
435 pm pitch leads to acceptable performance up to 
~175 GHz (~1.7 mm), which is the upper edge of our 
band. 

Comparing this relation between Vf and n eS to ana- 
lytic models [15,25] for the case where the electric 
field of the incident wave is perpendicular to the 
grooves, we found discrepancies that are reduced 
as the pitch decreases and we approach the quasi- 
static limit p/l«c 0.1. Given that our design does 
not operate in this limit and simple analytic models 
have insufficient accuracy for our purposes, we resort 
to numerical simulation to optimize the geometry of 
our AR layers. 

The variation in the electrical thickness compared 
to the physical thickness (Fig. 4, right) is shown for 
the case of a quarter-wave AR layer with index ^/n^. 
Additional simulations show that the electrical 
thickness depends on the index of refraction of the 
material on either side of the metamaterial dielectric 
layer. We account for this small effect in the final 
numerical optimization of the multilayer design. 

With these relations (Fig. 4) in hand we convert 
our two-layer coating into the desired parameters for 
the pillar geometry. In the final step we performed a 
numerical optimization of the AR layers at 15° angle 
of incidence using HFSS, with the constraint that the 
cut geometry must match the kerf geometry of com- 
mercially available dicing saws (see Section 6). This 
step ensures manufacturability, accounts for any 
discrepancies between the electrical and physical 
thickness implied by our numerical simulations, and 
improves the performance at larger angles of inci- 
dence. Table 1 gives the parameters for the resultant 
design, and Fig. 5 shows the simulated performance 
for this AR design as a function of frequency for 
several (0°, 15°, and 30°) angles of incidence. Even 
at 30° incidence the band averaged reflections are 
at -26 dB and differences between the two linear 
polarization states (TE and TM; see [ 34 ] for defini- 
tion) are below 0.5% for all frequencies. At 15° inci- 
dence, average reflections are below -31 dB. The 
tolerance of the design to manufacturing errors are 
described in detail in Appendix A. 

4. ACTPol Lenses 

The ACTPol lenses are cylindrically symmetric 
plano-convex designs in which the convex surface 


simulated performance 



modeled geometry 


200 urn 


Fig. 5. Left: Simulation of the performance of the ACTPol coating 
at a range of angles of incidence. Right: Unit cell geometry (based 
on measured cut profiles) modeled using HFSS. 


is aspheric — a conic section with four perturbing 
terms proportional to the fourth, sixth, eighth, and 
tenth orders of the distance from the axis. The design 
optimization procedure for the reimaging optics is 
similar to that described in [20], the primary differ- 
ence being that the ACTPol optics are required to be 
image-space telecentric to optimize the coupling to 
the planar feedhorn arrays. This was accomplished 
by constraining the chief rays at each field point to 
be near normal incidence at the focal plane and by 
allowing the tilts of the lenses and the focal plane 
to vary during the optimization. The resulting design 
achieves Strehl ratios greater than 0.93 across the 
150 GHz focal planes without accounting for the 
Gaussian illumination of the feedhoms, which effec- 
tively improves the image quality. All three of the 
ACTPol optics tubes use the same three silicon lens 
designs with the positions and tilts adjusted to opti- 
mize the coupling to the ACT Gregorian telescope [20]. 

The AR coating approach we have developed con- 
strains the perimeter of the lens designs to enable 
accurate clamping during both machining of the AR 
coatings and cryogenic cycling of the lenses. Figure 6 
shows a cross section of the perimeter of an ACTPol 
lens design. Each lens includes a handling ring along 
its perimeter that is not AR coated and is used for 
mounting the lens during machining as well as 
clamping to cool the lens in the cryogenic receiver. 
The outer corners of the perimeter are chamfered 
to minimize chipping. There are steps rising from 
the perimeter to the lens surface to provide clearance 
between the perimeter clamping region and the AR 
layers. The lens blanks were machined by Nu-Tek 


Table 1 . Parameters for the ACTPol Metamaterial Coating 


Parameter 

Symbol 

Dimension 

(pm) 

Dimension 
(units of p) 

Pitch 

P 

435 

1.0 

Kerf“ 


41 

0.094 


k 2 

190 

0.437 

Depth 

d\ 

217 

0.499 


d 2 

349 

0.802 

Volume fill 

v n 

0.820 


factor 




Vf2 

0.317 



“Nominal average values. 


ACTPol v28b final 1 50 GHz design - Lens 2 Edge 


8 mm 



Fig. 6. Cross section of the perimeter of the ACTPol lens 2, which 
is designed to enable machining of the AR layers and clamping the 
lens without damaging the AR coatings. 
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Precision Optical Corporation and achieved 5 pm 
tolerances. 

5. Properties of Silicon 

Silicon manufacturing can produce different grades 
of material, including ultrahigh purity silicon pro- 
duced by the float zone process [ 36 ] that has a neg- 
ligible density of impurities and silicon produced by 
the Czochralski process [ 36 ] that has a higher level of 
impurities. As charge carriers and associated states 
introduced by impurities are the cause of dielectric 
loss, use of the highest available purity (as inferred 
from the room temperature resistivity) silicon mini- 
mizes the dielectric losses. Measurements of the 
refractive index, dielectric permittivity, and loss 
tangent of various high purity and high resistivity 
silicon samples over a range of frequencies and tem- 
perature have been reported [37-39]. For example, 
ultrahigh purity silicon has been measured to have 
a loss tangent of ~1 x 10 -5 at room temperature 
which would correspond to <1% loss in the ACTPol 
optical system. Unfortunately, surface tension limits 
the zone-melt purification technique used to produce 
ultrahigh purity silicon to diameters below about 
200 mm. Therefore, the substrates available for large 
diameter lenses considered in this work must be fab- 
ricated from Czochralski silicon. For our low temper- 
ature application we expect the bulk conductivity 
of silicon to freeze out, dramatically reducing the 
dielectric losses; however, other loss mechanisms can 
persist. We validate this general picture of the dielec- 
tric loss in silicon by optically characterizing the 
influence of the bulk resistivity as a function of fre- 
quency and sample temperature. 

Samples produced by the Czochralski process 
and readily available in 450 mm diameter stock with 
bulk resistivities specified to be in the range of 1 to 
>500 G-cm were characterized both at room temper- 
ature and 4 K using a Bruker 125 high-resolution 
Fourier transform spectrometer (FTS) with a Oxford 
Cryostat CF continuous liquid helium flow sample 
chamber. The cryostat is equipped with pair of 75 pm 
thick polypropylene windows that enable spectral 
measurement while allowing the sample to be held 
at a regulated temperature. The samples were cut 
to have a typical thickness of 180 pm, double side pol- 
ished, and placed in a 25 mm diameter optical test 
fixture at the focus of an f / 6 beam. The reflections 
from the two surfaces of the sample form a Fabry— 
Perot resonator for which the modeling is relatively 
simple, permitting measurement of dielectric proper- 
ties. The silicon samples are boron doped (p-type) to 
adjust the resistivity. The >500 G-cm resistivity sil- 
icon used in the ACTPol lenses has the minimum 
dopant level. 

Each sample’s transmission was measured between 
240 GHz (8 cm -1 ) and 18 THz (600 cm -1 ) using differ- 
ent combinations of sources, beam splitters, and 
detectors for three frequency bands of 240-450 GHz 
(8-15 cm -1 ), 450-2850 GHz (15-95 cm -1 ), and 
2.85-18 THz (95-600 cm- 1 ). 


The spectral resolution employed, 7.5 GHz 
(0.25 cm -1 ), fully resolves the sample’s spectral fea- 
tures. Sliding stages permit the sample or a reference 
clear aperture to be moved into the FTS beam while in 
the cryostat for in situ calibration. 

The transmission spectrum was modeled as a series 
of homogeneous plane parallel dielectric layers [ 40 ] . 
The dielectric function for the silicon was approxi- 
mated by a classical Drude dispersion model [41]: 


£*(<u) = £ 


* 

OO 


&>n 


co * (co T ir) 


( 2 ) 


where = e' r + ie" is a complex function of frequency 
co, the damping rate T, and the contribution to the 
relative permittivity of higher energy transitions. 
The plasma frequency and the damping rates are re- 
lated by ojp = T/ e Q p, where e 0 is the permittivity of 
free space and p is the sample’s bulk resistivity [42] . 
In this approximation, the material’s free carriers 
are treated as classical point charges undergoing 
random collisions, and the resulting damping is as- 
sumed to be independent of the carrier energy. We find 
this representation suitable to represent the sample’s 
properties over the spectral and temperature ranges 
of interest. 

The sample thickness is known to ±0.5 pm at room 
temperature and corrected for thermal contraction 
[ 43 ] . Since the fringe rate is proportional to the prod- 
uct of the refractive index and sample thickness, un- 
certainties in the sample thickness directly limit the 
precision of the determination of (ftfe^). Measure- 
ments of an optically polished crystalline quartz 
sample with accurately known thickness were used 
to measure a 1% amplitude uncertainty across the 
entire FTS band. This calibration uncertainty leads 
to a corresponding reduction in the measurement’s 
sensitivity to (sfe^). These uncertainties were ac- 
counted for in fitting the FTS data. These fits pro- 
duced root-mean-square deviations in the range of 
~0.005 to 0.016 between the model and the observa- 
tion spectra. Representative data for the p > 500 and 
1 < p < 5 G-cm samples are shown in Fig. 7 for sam- 
ple physical temperatures of T pliys = 297 and 4.5 K. 

The results of these fits are shown in Table 2. As 
anticipated the finite plasma frequency and damping 
rate contribute significantly to the room temperature 
losses; however, the p > 500 G-cm sample could be 
approximated by a dielectric constant over the spec- 
tral range of interest (e.g., T and co p are consistent 
with zero). The influence of free carrier collisions at 
room temperature on the optical response becomes 
more pronounced in the 1 < p < 5 G-cm sample, 
manifesting as a reduction in transmission at low 
frequencies. However, as the 1 < p < 5 G-cm sample 
is cooled, bulk conduction is suppressed, and this 
effect shifts to lower frequencies than those of inter- 
est for millimeter-wave applications. The observed 
dielectric parameters are a weak function of temper- 
ature below ~30 K, as anticipated given the exponen- 
tial thermal dependence of the bulk resistivity [ 44 ] . 
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Silicon: 1 < p < 5 (£2-cm) Silicon: p > 500 (£2-cm) 



024602468 


Frequency (THz) Frequency (TFIz) 

Fig. 7. Measured silicon transmission at room temperature 
(upper panels) and 4.5 K (lower panels). The left and right columns 
in the panel are for samples specified as 1 < p < 5 Q-cm and 
p > 500 fi-cm, respectively; the sample thickness is ~180 pm. 
The figure contains the measured FTS data (red), model (dashed 
blue), and residual (green). 


Both samples show ~1% shifts in the magnitude of 
the index of refraction in « -y/eQ upon cooling. This 
shift is accounted for in our optical design. At low 
temperature both the high and low resistivity sample 
give an upper limit for the loss tangent of 
tan 5 ~ e'^/e'oo < 7 x 10" 5 . For the ACTPol optics at 
150 GHz, this corresponds to an absorptive loss of 
<5% for our optical system comprising three lenses. 
The ACTPol lenses were fabricated from p > 
500 G-cm silicon as it may reduce loss and facilitates 
room temperature testing. Our models and observa- 
tions suggest that the carriers will be frozen out at 
4 K for both high and low resistivity silicon, but 
given the small cost differential high resistivity 
silicon is the natural choice. We note that the con- 
ductivity of silicon at room temperature, and there- 
fore the loss, is a strong function of UV irradiance. 
We have not investigated this dependence cryo- 
genically. 


repeatedly cut with micrometer-level precision and 
cutting speeds up to several centimeters per second. 
Unfortunately, commercial silicon dicing machines are 
not designed for curved surfaces and cannot accommo- 
date the large diameter lenses required for this 
project. Therefore we constructed a custom three-axis 
silicon dicing system to fabricate these structures. 

The fabrication system is shown in Fig. 8. It con- 
sists of a micrometer accurate three axis stage on 
which we mount an air bearing dicing spindle and a 
micrometer accurate depth gauge. The spindle and 
retractable depth gauge are attached to the vertical 
stage (z axis), which rides on the horizontal (y axis) 
stage. A lens is mounted on an aluminum mounting 
plate on a horizontal stage (x axis) below the spindle. 
This plate permits the lens to be accurately rotated 
by 0°, 90°, 180°, and 270° and carries a reference wa- 
fer that is used in setting up the blades prior to cuts. 
A side-looking microscope is mounted parallel to the y 
axis to characterize test cuts on the reference wafer. 
Flood cooling water is sprayed on the dicing blade 
while cuts are made to carry away debris. A temper- 
ature controlled water bath is used to regulate the tem- 
perature of the spindle and maintain the flood coolant 
and air surrounding the dicing system within 1°C. 

Our system does not show any appreciable change 
in cut shape or surface damage when cutting at the 
maximum travel speed of 50 mm/s. We conserva- 
tively operate at 25 mm/s, for which it takes a total 
of 12 h of machining per lens side. Factoring in setup 
time it takes six 8-h days to fabricate a single lens. 

Figure 9 shows photographs of one of the ACTPol 
AR coated silicon lenses. The fabricated features are 
sufficiently robust to permit handling the lenses by 
touching the AR coatings. The manufacture resulted 
in fewer than 10 out of 500,000 posts with damage 
visible by eye. The shape of the cut profiles was 
evaluated by making cuts on a reference silicon wafer 
before cutting (precut) and after cutting (postcut) the 
lens and measuring them using the side-looking mi- 
croscope. In addition, reference cuts were examined 
at frequent intervals during the process of cutting 


6. Fabrication 

The AR coatings were fabricated using silicon dicing 
saw blades, which are available in widths ranging 
from 10 to hundreds of micrometers. Mounted on a 
commercial silicon dicing saw machine, these blades 


Table 2. Silicon Drude Model Fit Parameters 
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[Q-cm] 
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4 
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[THz] 

m p /2ir 

[THz] 

>500 

297 

11.7(1) 
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— 

— 


4.5 

11.5(5) 

<0.0008 

— 

— 

1 to 5 
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0.0046 
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11.5(52) 

0.0049 
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100 
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0.0028 
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Fig. 8. The custom, micrometer accurate, three axis silicon dicing 
system used to fabricate AR coatings on the lenses. The labels 
identify the key components and axes described in Section 6. 
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Fig. 9. Left: Photograph of the curved front surface of an AR coated ACTPol silicon lens (lens 2) with zoomed view of a small patch. Right: 
Photographs taken with side-looking metrology microscope of a reference wafer taken prior to and after cutting the coatings on a lens. 


the lens, and the whole range of profiles were found 
to be within acceptable tolerance limits. Comparison 
of the precut and postcut measurements (see Fig. 9) 
show that the wide blades cut repeatably with 
negligible evolution to cut profile while the narrow 
blades show some evolution in width in the upper 
third of the cut. Simulations show that this evolution 
leads to a few tenths of a percent increase in the 
reflectance. The accuracy of the depths of the cuts is 
limited by a 3 pm uncertainty in the zero point for the 
blade depth and a ±7 pm surface error in the lens. 
This surface error is due to imperfections in the 
lens mounting plate that stress and warp the lens. 

The reflectance of the completed ACTPol lenses 
were measured using a scalar reflectometer. This 
reflectometer consisted of a tunable narrowband 
continuous wave source which illuminated the flat 
side of lens through a 20° full width half power horn, 
from a distance of 5 cm, and tilted at an angle of 15° 
relative to the lens surface normal. A receiver con- 
sisting of an identical horn coupled to a detector 
diode was placed at the mirror image of the source 
horn relative to the plane defined by the point of 
maximum illumination of the lens and perpendicular 
to a line joining the source and receiver. This system 
was calibrated by (1) placing an aluminum reflector 
at the same position as the lens to normalize the 
peak reflection to unity, and (2) by removing the lens 
and calibration reflector to measure the stray reflec- 
tions which were found to be negligible. 

The results of this reflectance measurement are 
shown in Fig. 10 for both the TE and TM polariza- 
tions. Precise modeling of this measurement configu- 
ration requires accounting for the interference 
between reflections from the flat and curved lens 
surfaces, which is beyond the scope of this work. 
However, the results are expected to be intermediate 
between the reflectance from a single AR coated 
surface (single sided) and the easily simulated result 
for two flat coated surfaces separated by the central 
thickness of the lens (two sided). These two cases are 
presented in Fig. 10 for both linear polarization 
states. These simulations incorporate the measured 


Lens Reflectance at 1 5° incidence 



120 130 140 150 160 170 

Frequency (GHz) 


Fig. 10. Comparison between simulations and measurements 
of the fabricated lens. The two-sided simulations model the AR 
coatings on both sides of a flat silicon sample of thickness equal 
to that of the lens at its center. The single-sided simulations model 
the AR coatings on only one side of a flat silicon sample. 

cut profiles shown in Fig. 9. We empirically found 
that moving the lens so that it is well centered on 
the beam increases the interference effects (TM case) 
while moving the lens off center reduces these effects 
(TE case). Based on the reasonable agreement be- 
tween these simplistic simulations and the measure- 
ment we are confident we are reducing reflections to 
a few tenths of a percent over the range of angles of 
incidence required for ACTPol. 

7. Conclusion 

We have described a new approach for AR coating 
silicon lenses over broad bandwidths and a range 
of angles of incidence. Simulations, backed up by 
measurements of an AR coated lens, show that the 
fabricated coatings presented here reduce reflections 
below a few tenths of a percent between 125 and 
165 GHz for angles of incidence between 0 and 30° 
for cryogenic applications. We have developed a 
micrometer accurate three axis silicon dicing saw 
facility and are using it to manufacture AR coated 
lenses for ACTPol. We have shown that a range of 
p-type silicon doping levels can achieve low loss at 
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cryogenic temperatures using silicon samples from 
boules as large as 45 cm diameter. This approach for 
implementing wide-bandwidth AR coated silicon 
lenses is applicable for millimeter and submillimeter 
wavelength ranges and can be expanded to wider 
bandwidth by adding additional layers. 

Appendix A: Tolerances 

The manufacturing tolerances for the eight types of 
machining errors shown in Fig. 11 were evaluated 
based on a fiducial design. These include (1) errors 
in the pitch, Sp; (2) errors in the centering of the 
deeper groove relative to the shallow groove, Sc; (3) 
errors in the kerf widths, Sk 1 , Sk 2 ; (4) slopes in the 
walls of the grooves parameterized by <5si and $S2> 
which represent the distance the upper and lower 
grooves tilt inward at the groove bottom if the mean 
width is fixed; (5) errors in the depth of either groove, 
Sdi, Sd 2 ; (6) differences in the depth of grooves in 
the two orthogonal directions assuming the correct 
mean depths are parameterized <5o? LRl for the upper 
and <5 cZ L r 2 for the lower groove; (7) fabricating layers 
on a surface tilted at an angle <p (see Fig. 11) and 
where the depth of each groove is measured at the 
center of the groove (e.g., the effect of applying these 
coatings at the edge of a curved lens); and (8) rows 
of broken posts. The tolerance to variation in the re- 
fractive index of the silicon substrate, <5nso was also 
evaluated. The tolerances for these parameters were 
quantified based on a fiducial model similar to what 
was fabricated but with straight-walled pillars using 
HFSS simulations in which the parameters in each 
family were varied separately. For the majority of 
these effects the results were distilled to a band 
averaged reflection at 15° angles of incidence. The 
sensitivity to each parameter was quantified as 
the displacement needed to bring about a 3 dB 


types of machining errors 

pitch centering kerf sloped walls 


1 j' 1 j 

n n 


f ] 

L J - 1 j 

L 1 r J L 1 oi i*. 


L 1 ■ L Hh- 

j ! ! 

! ! ; ; s 

k 2 /2 

\ \ 5s 2 

Li Li 

8p 

8c 8k! 

12 

-HI— 

8si 


depth asymmetric depth tilted surface damaged row 




Fig. 11. This figure shows the families of machining errors con- 
sidered in the tolerance analysis and described in the text. The 
dashed lines show the target shape for the grooves and that for 
the pillars while the white space in the gray represents the 
material as actually cut. Therefore difference between the dashed 
region and the white region represent machining errors. For the 
asymmetric depth and damaged row errors the sketches show 
the geometry of the errors. 


increase in reflectance. For the case of tilted surface 
and for broken posts the results of a small number of 
simulations were evaluated and compared to the fi- 
ducial performance. In the remainder of this section 
we discuss these results. 

A. Pitch 

Variations in the pitch change the effective index of 
both layers. Figure 12 shows the band averaged 
reflectance as function of variation in the pitch. A 
conservative tolerance of <5p < 3 \im(Sp/p < 0.0075) 
limits the change in reflectance to be less than 3 dB. 

B. Centering 

The performance was insensitive to errors in the 
centering of the deeper grooves within the shallower 
grooves. This is consistent with both layers of pillars 
behaving as layers with a constant effective index 
of refraction. Aperiodic misalignments of the deeper 
grooves within the shallow grooves are equivalent 
to variation in the pitch. 

C. Kerf 

The upper left panel of Fig. IB presents the impact of 
errors in the kerf. Controlling Ski an d Sk 2 to 3 pm en- 
sures a less than 3 dB degradation in performance. 
Varying both parameters in the same direction 
(e.g., Ski — has little impact on the performance. 
We have chosen the fiducial design to be near the 
high end of the most favorable region to make the 
manufacture relatively immune to blade wear, which 
could have (but did not) narrow the kerf width as 
machining proceeds. 

D. Sloped Walls 

The upper-right panel of Fig. IB presents the impact 
of slopes in the walls left behind by the dicing saw. 
Conservative estimates for the 3 dB tolerance for 
these two parameters are 6 pm for <5si and 3 pm for 
Ss 2 . Sloping both layers by a similar amount results 
in negligible degradation. 


variation in pitch 



Sp (um) 


Fig. 12. Band averaged reflectance at a 15° angle of incidence as 
a function of changes in the pitch. The gray region shows range for 
which the performance is within 3 dB of the fiducial design. 


20 December 2013 / Vol. 52, No. 36 / APPLIED OPTICS 


8755 


variation of substrate index 




6d, [nm] 


-30 -15 0 15 30 

6d LRi [pm] 


-6-3 0 3 

5k, l|jm] 

depth error 


0 3 6 9 12 

5s, [pm] 

asymmetric depth 


kerf error 


sloped walls 


cu ~a 
c 2 

16.0 - £ 
•— <u 
14.5 | | 
<T3 E 
13.0 -c o 


1.5 

0.0 


Fig. 13. This figure shows the impact of errors in the kerf, slopes 
in the sidewalls, overall depth errors, and asymmetric differences 
in depth between the two orthogonal cuts. The color scale presents 
the degradation in reflection at 15° incidence compared to the sim- 
ulations in Fig. 5. The fiducial design is at the coordinate (0, 0) in 
all the plots as highlighted by a “+” or arrow. The outer edge of the 
light blue contour (see label) represents 3 dB degradation. The 
horizontal (vertical) axis represents errors in the inner (outer) 
layer for the labeled parameters. 


E. Depth 

The lower left panel of Fig. 13 presents the impact of 
depth errors. The 3 dB tolerances for these parame- 
ters are 7.5 pm for 8di and 10 pm for 5d 2 - Expressed 
in terms of optical path length the sensitivities are 
identical. 

F. Asymmetric Depth 

The lower right panel of Fig. 13 presents the impact 
of asymmetric depth. The 3 dB tolerances for these 
parameters are 15 pm for <5 g? LRi and 12 pm for &7 L r 2 
for asymmetry in directions common to both grooves. 
This error also affects the 90° rotational symmetry 
of the geometry resulting in birefringence, but the 
cross polarization can be no larger than the absolute 
reflection. 

G. Application to Tilted Surfaces 

The left panel of Fig. 14 compares the performance of 
these layers cut on a surface tilted by 10° compared 
with the performance when cut on a flat surface. 
Fortuitously, this effect improves performance at 
large angles of incidence and has little effect at small 
angles. Applying this coating by making grooves par- 
allel to the lens symmetry axis of a curved lens will 
produce acceptable performance. 

H. Groups of Broken Posts 

Assessing the impact of the broken posts is not 
straightforward because the region and extent of 
illumination varies between the three lenses. The 
degradation due to broken posts would thus be a 
function of the effective aperture area of the lenses 


CD 
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surface tilted by 10° 
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Fig. 14. Left: Simulated impact of having the surface of the lens 
tilted by 10°. Right: Simulated impact of variation in the refractive 
index of the silicon substrate. 


and the location of the defects in addition to the total 
area affected. Accurately quantifying this is beyond 
the scope of this work. Qualitatively, the degradation 
would be negligible if the fraction of total area 
affected is small and, more importantly, the length 
scale of the individual affected areas is smaller than 
the wavelength of incident light. Simulating the 
impact of broken posts is computationally difficult. 
Therefore we resorted to a simplified model consist- 
ing of a unit cell containing seven posts in a row 
with one broken off completely. This is equivalent 
to having one-seventh of the posts destroyed. This 
pessimistic case produces band averaged reflections 
below 2%, which is acceptable, though it does produce 
nontrivial cross polarization. We used simulation to 
estimate that keeping the number of broken posts 
below 1/700 would result in a negligible degradation 
in the overall performance of the AR coating. 

I. Index 

Figure 14, right, shows the band averaged reflectance 
as a function of variation in the index. A conservative 
tolerance for the index is <5re S i < 0.1 (Snst/ngj < 0.03). 
This leads to increases in the reflectance of less than 
3 dB over the fiducial design. This tolerance is weaker 
than that imposed by the optical design. 

The 3 dB sensitivities derived in this section are 
summarized in Table 3. The band averaged reflectiv- 
ity of the fiducial design is -31 dB at 15° incidence, 


Table 3. Summary of the 3 dB Sensitivities for All Parameters 
Described in the Text 


Machining Error 

Parameter 

Sensitivity 

(pm) 

Tolerance 
(units ofp) 

Pitch (p) 

Sp 

3 

7.5 x 10 -3 

Centering 

Sc 

See text 


Kerf 

8k i 

3 

7.5 x 10 -3 


8k 2 

3 

7.5 x 10“ 3 

Sloped walls 

(5s i 


1.5 x 10“ 2 


Ss 2 

3 

7.5 x 10 -3 

Depth 

Sd 1 

7.5 

1.9 x 10“ 2 


Sd 2 

10 

2.5 x 10“ 2 

Asymmetric 

drfi.R, 

15 

3.8 x 10“ 2 

depth 

(5c?lr 2 

12 

3 x 10“ 2 

Tilted surface 

<P 

See text 


Damaged row 

r 

-1/700 

NA 

Index 

Sn 

0.1 

NA 
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which is significantly better than the requirements 
for the system. 
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